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Abstract 

The State-Specific Kohn-Sham Density Functional Theory arXiv:physics/0506037 is used to 
derive the Kohn-Sham exchange-correlation potential v xc and exchange-correlation functional E xc 
as explicit functionals of v s and (fx, where v s is the local, one-body potential from the Kohn-Sham 
equations, and <f>\ is the spinless one-particle density matrix from the Kohn-Sham noninter act- 
ing state, say \<pi)- In other words, is the ground state eigenfunction of the noninteracting 
Schrodinger equation with the one-body potential v s . For simplicity, we only consider noninter- 
acting states that are closed-shell states and interacting states that are nondegenerate, singlet 
ground-states. 
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I. INTRODUCTION 



The Kohn-Sham version of density functional theory plays a major role in both quantum 
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chemistry and condensed matter physics [1|, |2j, |3|, y, |5|, |fj |7| . Unfortunately, the exchange- 
correlation functional E xc is an unknown, implicit functional, and there is no systematic 
method to improve approximations. Recently we have derived a generalization of the Kohn- 
Sham approach in which the correlation energy E co is assumed to be an explicit functional 
of v and p±, where v is the external potential from the interacting target-state, and p\ is the 
spinless one-particle density matrix from the noninteracting states 0,0]- In this approach, 
errors from Coulomb self-interactions do not occur, nor the need to introduce functionals 
defined by a constraint search. Furthermore, the exchange energy E x is treated as in Hartree- 
Fock theory, as an explicit functional of p\. Below, we use this approach to derive the Kohn- 
Sham exchange-correlation potential i> xc and exchange-correlation functional E xc as explicit 
functionals of v s and cpx, where v s is the local, one-body potential from the Kohn-Sham 
equations, and (fx is the one-particle density matrix from the Kohn-Sham noninteracting 
state, say \(fix)- In other words, \<pi) is the ground state eigenfunction of the noninteracting 
Schrodinger equation with the one-body potential v s . 

II. STATE-SPECIFIC KOHN-SHAM DENSITY FUNCTIONAL THEORY 

In state-specific Kohn-Sham density functional theory 0, 0] , we use an energy functional 
E v [pi] that is assumed to be an explicit functional of the external potential v and the spinless 
one-particle density matrix p±, where p\ comes from a closed-shell determinantal state, say 
this energy functional is given by 

/vj> \H 1^ ) 

rp r i _ \*vp 1 \ 11 v\*vp 1 / /-.x 

M ~ (V I* ) ' () 

where the trial wave function ^ vpi generates the exact, or target, wave function, say ^/ n , 
under the following conditions: 

|*„ w ) = |* n ); Qi — n — >N,v, (2) 

where n is the electron density of \I/ n , and is the ground state singlet eigenfunction of 
H v . The right side notation of Eq. indicates that £i is a spin- less one particle density 
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matrix that delivers the density n; according to the Hohenberg-Kohn theorem 
n also determines the external potential v , and n determines the number of electrons N. 
Using Eqs. JTJ and (J2J), we have 

E v [gi]=S n , qi — >n, n — ► N,v. (3) 

Here, £ n is the exact electronic energy of the target state: 

H v \^ n ) = £ n \q n ); n-^N,v. (4) 

where the Hamiltonian operator is given by 

H v = f + V ee + V v , (5) 

and we have 

T = D-|V?), (6) 



i 

The electronic energy functional can also be written as 

£,M = J *i [-I v ?/ 9 i( r i' r 2)] r2=ri + y drv(r)p a (r) 

+ £j[p s ] + £ x [ Pl ] + £ co [pi,v] + / <2™(r)p c (r), (9) 
where the Coulomb and exchange energies are given by the following: 

Ej[Ps] = \J J rr 2 1 rfr 1 rfr 2 p(r 1 )p(r 2 ), (10) 



-E x [pi) = J r 12 1 rfrirfr2pi(ri,r 2 )pi(r2,ri), (11) 
and the correlation-energy functional is 

^co[p^] = - - - (pi|T| Pl ) + (pil^lPi). (12) 

Furthermore, p c is the correlation density of the trial wave function, i.e, we have 

p c (r) = <y( r )i*«pi> _ Ps(r) = a _ Ps(r); ^ _> fij P1 _ Ps; (13) 



and f is the density operator; n is the density of *& vpi ; p s is the density of p\. Since p c ( r ) is 
a functional of v and pi, we can also write p c [pi,v](r). 

A determinant al state with the density matrix g 1 satisfies the following noninteracting 
Schrodinger equation: 

f^^ 1 (r,)|^ 1 ) = 2^ £ J \ Ql ), (14) 

1=1 \ w / 

where the generalized, or exact, Fock operator T ei is given by 

K = -\V 2 + v + v n J + Of + t>£ + . (15) 
Here, the Coulomb operator is defined by 

(ri)x(ri) = J dr 2 rr 2 V(r 2 )x(ri), (16) 



and the exchange operator Of, correlation operator 0£J and external-correlation operator 
are defined by their kernels: 

pu \ SE x [p u v] 1 _! 

<(r " r2> = ^M = -2 r " ft(r " r,) ' (17) 

- <(/ £y ' <»> 

Our energy functionals -E„ are implicit functionals of the noninteracting density p s . Hence, 
any one-particle density-matrix that yields the interacting density minimizes our energy 
functional, i.e., we have 

S n = E v [ Ql \ = E v [q[] = E v [q'[\ ■■■ , (20) 

where 

ra(r) = gi(r, r) = ^(r, r) = ^'(r, r) • • • . (21) 

Assuming n is a noninteracting f-representable density, there exist a noninteracting state, 
say \<pi), that has n as its density: 

n(r) = ^(r,r), (22) 



and this determinant — assuming it is a closed-shell determinant — is the ground-state solution 
of the following noninteracting Schrodinger equation: 

£/( r «)l^> = 2 (£ e ») M, (23) 

i=l \ w J 

where 

f = -\V 2 + v S) (24) 

and v s is a local potential. Therefore, the canonical occupied orbitals from \ipi) satisfy the 
following one-particle Schrodinger equation: 

f<p w = (~|V 2 + v + Vj + f xc ) <p w = e w cj) w , <j) w e<pi, (25) 

where, with no loss of generality, we have required v s to be defined by 

v s = v + v 1 } + v xc . (26) 



Using the approach by Sala and Gorling 



111 ] , but permitting the orbitals to be complex, 



it is readily demonstrated that t> xc is given by 

^xc(r) = 2~^j J rfr i [2^(ri,r)^!(r,ri) - (/^(r, r x ) J dr 2 <Pi(r 2 , t)w(ti, r 2 ) (27) 

+V5i(ri,r)( / 9 1 (r,r 1 )v xc (ri)] , 

where to is the kernel of the nonlocal potential w pi , given by 

+ + (28) 

and these operators appear in the exact Fock operator given by Eq. (|15|) . By substituting 
v xc repeatedly on the right side, we can obtain an expansion for v xc : 

^xc(r) = ^yjl 2w ( r ^ r Vi( r , r i) - ¥>i(r, ri)^i(r 2 , r)w(r 1; r 2 ) 

+ v?i(ri, r)v?i(r, n)— ; ^-{u;(r 2 , r^tp^rx, r 2 ) - -y>i(ri, r 2 )^ 1 (r 3 , r x )-u;(r 2 , r 3 )} 
n(rij 2 

+ y>i(ri,r)y?i(r,ri) - 1 y?i(r 2 , ri)y?i(ri, r 2 )— j-r-w(r 3 , r 2 )y?i(r 2 , r 3 ) H ], (29) 

ln{Yi) n{r 2 ) 

where in this relation it is understood that there are integrations over the dummy variables 
r 1; r 2 and r 3 . The leading term of Eq. (}2"9"j) is the Slater potential [3, [3, Q|; this term 
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also appears within the Krieger-Li-Iafrate (KLI) approximation of the optimized potential 
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The orbitals <p w satisfying Eq. (|25|) are the Kohn-Sham orbitals |l8j : |v?i) is the Kohn- 
Sham noninteracting state. However, / differs from the Kohn-Sham operator, since, in 
addition to depending explicitly on tpi, instead of n, f depends explicitly on the external 
potential v from the interacting Hamiltonian H v . Furthermore, the external-correlation 
operator {ffi does not appear in Kohn-Sham formalism. And, unlike the original Kohn-Sham 
approach [181 ] . the iV-represent abili ty p roblem does not arise, nor the need to introduce a 

22| to avoid this problem. 



constraint-search definition [13, |20j, 
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III. THE KOHN SHAM EXCHANGE CORRELATION POTENTIAL 

According to Eqs. ([THjl. and (fTUjl. w pi is a functional of p\ and f , indicating that we 
can, symbolically speaking, represent Eq. (j2HJ) by 

v xc (r) =v xa [(p 1 ,v](r), ifx — >n — >v. (30) 

In other words, v xc is also a functional of v and ip\, where ip\ determines n, and n deter- 
mines v. 

Note that (pi and v from Eq. (J3U|) are not independent. Since, for a given v, the one- 
particle density matrix ipi which determines v xc from Eq. ([$()[). and then gives f s from 
()26[) . must also be the one-particle density matrix ipi from the noninteracting state \ipi) 
that satisfies Eq. ()23|1. However, for a given or f s , we can also think of v as a dependent 
variable to be determined. In that case, we choose v s , construct / using Eq. (|24p . and obtain 
the one-particle density matrix, say ipi, that satisfies Eq. (j23|) . The external potential v is 
then a simultaneous solution of Eqs. (|2*B|) and ([3L)p. 

Substituting Eq. (JSDJ) into Eq. JJSJ, gives 

u = u s - - v xc [cpi,v], N,v s — > <px — >n — > N,v,v s , (31) 

where the notation on the right side indicates that iV and v s determine ipi, as indicated by 
Eqs. d2ni and (JUJ); also, ip\ determines n; n determines N, v, and v s . 

By substituting v on the right side of Eq. ([3~T|) repeatedly, as in Eq. (j29jl . we can remove 
it, giving, symbolically speaking, 

u = t#i,uj, (32) 
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Hence, this relation gives the external potential v from the interacting system ty n , for n — > 
v, as a functional of the local potential v s and the one-particle density matrix tpi from the 
noninteracting state \ipi), where \tpi) is an eigenfunction of the noninteracting Hamiltonian 
given by Eq. (J2HJ), and this noninteracting Hamiltonian is defined by the local potential 
v s , as indicated by Eq. (|24|) . Furthermore, the noninteracting state \tpi) shares the same 
density with \I/ n ; \(fi) is the Kohn-Sham determinantal state when n is noninteracting v- 
representable. 

Since there is a one-to-one correspondence between local potentials v s and noninteracting 
ground states Eq. (f3~2|) can be written as 



v = v[(p u v a ], v s < — > (fx, 



(33) 



where the right side indicates the one-to-one correspondence. 
Substituting Eq. (J33j) into Eq. (|H1) gives 

^xc = V xc [(fl , V s ] , V s < > If i . 

Denoting v^ the Kohn-Sham exchange- correlation potential [l| 
for non-interacting i>-representable densities, we have 



(34) 
1811 . at least 



^xcfM = v xc [(fi,v s ), n — ><pi,v 8 , v s < — > <fi. 



(35) 



In order to obtain the density functional on the left side of Eq. (1351) , we need the function- 
al fi[n] and v s [n). According the to Hohenberg-Kohn theorem l|, l2,[l0(, for noninteracting 
v-representable densities, these functionals exist. For other densities, however, these expres- 
sions must be generalized by, for example, using some modified constraint search definition, 
or, perhaps, by an approach that permits v s to be nonlocal. 



IV. THE KOHN-SHAM CORRELATION ENERGY FUNCTIONAL AND EX- 
CHANGE CORRELATION FUNCTIONAL 



Using t 
theorem 



Li 



re notation form Eq. (jljl . the universal functional from the Hohenberg-Kohn 
, can be written as 



F[n] = (* n \[f + V ee )\* n )i n^N,v. 



(36) 
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where \l/ n is a singlet, ground state wave function. Previously we have shown that the 



23, 



2d, Urn can be written as an 



correlation energy from many body perturbation theory 
explicit functional of v and p\ [26(. In a similar manner, but using less restrictive energy 
denominators, the universal functionals F can be shown to be an explicit functional of v 
and pi js?!, where this functional does not implicitly depend on pi, i.e., any that has 
considerable overlap with |\& n ) can be used as a reference state. So, we can write 

F[n) — F[v,px], n^v, n, Pl —>N. (37) 

For noninteracting v-representable densities n, we can use Eq. (f3*2~j). giving 

F[n] = F[(pi,v s ,pi), n — > <fii,v s , v s < — > <px, p\ — ► N. (38) 

Setting pi = cpi, we get 

F[n) = F[<pi,v„], n — ><pi,v„, v s < — > ip 1 . (39) 

Hence, assuming the existence of the explicit functional, given by Eq. (|37j). Eq. (|3*9*jl gives 
this functional as an explicit functional of v s and (fx, where v s is the local, one-body potential 
from the Kohn-Sham equations, and (fx is the spinless one-particle density matrix from the 
Kohn-Sham noninteracting state \ fx)- 

Using Eq. ()39)1 . the exchange-correlation functional from the Kohn-Sham formalism 

is 

E™[n} = F[<px,v s }-Ej[n)-{<px\f\ipx), n^cpx,v s . (40) 
It is well know that the Kohn-Sham exchange functional is given by 0, 

£ x KS [n] = EM, n^yx. (41) 
Hence, the Kohn-Sham correlation-energy functional is 

E co [n] KS = F[ipx,v s ] - (<px\T\ipx) ~ (<pi\V<*\<Pi) , (42) 

where 

(<Pi\Vce\<Pi) =Ej[n] +E x [( Pl ]. (43) 
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APPENDIX A: THE KOHN SHAM CORRELATION-ENERGY FUNCTIONAL 
FROM STATE SPECIFIC KOHN SHAM DENSITY FUNCTIONAL THEORY 



Consider the energy functional E v [pi,v'} that is obtained by permitting the external 
potential from the trial wave function, say v', to differ from the one from the Hamiltonian 
H v ; generalizing Eq. (0) in this way, we have 



and Eqs. (0), ©, and Q become 

\^v> ei ) = |*n); v' = v, 6l — >n, n — ► N,v, (A2) 

K[Qi,v] = S n . (A3) 

K[pi, v'] = £i[pi, v) + E co [pi,v r \ + J drv(r)p c [pi,v')(r), (A4) 

where we have explicitly mentioned the dependence of p c on p\ and v', and the energy 
through the first order, Si, is given by 

Si[pi,v] = (pi\H v \ P i) = J dn [-IV? Pl ( ri ,r 2 )] r2=ri (A5) 
+ J dri u(ri)p(ri) + ^ J J dr 1 dr 2 r 12 1 p(ri) p(r 2 ) - J J J dvidr 2 r^ pi{r x , r 2 ) pi(r 2 , r x ). 

Using Eq. (J32J), we have 

E v [p 1 ,v'[<p , v v s ]]=S 1 [p l ,v]+E c0 [p 1 ,v'[i P ' 1 ,v s ]] + J drv^p^v'^v'Mr). (A6) 

Setting pi = <p' x 

E v Wi , i/„] = £j [ip'i , u] + £ co Wi , <] , ( A7) 
and suppressing the primes, we have 

E v [cpi , v s ] = Si [cpi , v] + E co [tpt ,v s ], ( A8) 

where we used 

p c [<p[,v'[<p[,v' s }} = 0. (A9) 
9 



Comparing Eq. ()A8|) with the energy functional from the Kohn-Sham formalism, for non- 
interacting f-representable densities, we have 

E* s [n] = E v [ip x ,v a \, (A10) 
Ef*[n] = E m [ip u v 8 ], (All) 

where n — ► <pi,v s , v s < — ► cpi. As in v^, in order to obtain the above density functionals, 
we need <pi[n] and yJn]. Recall that it is well know that the Kohn-Sham exchange functional 
satisfies Eq. (JUj) 
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